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We prove that, in the large-dimension limit, the high-density correlation energy Ec of two opposite- 
spin electrons confined in a _D-dimensional space and interacting via a Coulomb potential is given by 
Ec ~ —1/{SD^) for any radial confining potential V{r). This result explains the observed similarity 
of Ec in a variety of two-electron systems in three-dimensional space. 
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Understanding and calculating the electronic correla- 
tion energy is one of the most important and difficult 
problems in molecular physics. In this pursuit, the study 
of high-density correlation energy using perturbation the- 
ory has been particularly profitable, shedding light on the 
physically relevant density regime and providing exact re- 
sults for key systems, such as the uniform electron gas \V 
and two-electron systems [2]. The former is the corner- 
stone of the most popular density functional paradigm 
(the local density approximation) in solid-state physics 
[3]; the latter provide important test cases in the devel- 
opment of new explicitly correlated methods [H |S] for 
electronic structure calculations [F. Atomic units are 
used throughout. 

The high-density correlation energy of the helium-like 
ions is obtained by expanding both the exact [7] and 
Hartree-Fock (HF) [8 energies as series in 1/Z, yielding 



E{Z,D,V) ^E^^\D,V)Z^ + E^^\D,V)Z 

EnF{z, D, V) = £;(") {D, v)z^ + £;(!) {d, v)z 
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where Z is the nuclear charge, D is the dimension of the 
space and V is the external Coulomb potential. Equa- 
tions (IT]) and ([2]) share the same zeroth- and first-order 
energies because the exact and the HF treatment have 
the same zeroth-order Hamiltonian. Thus, in the high- 
density (large-Z) limit, the correlation energy is 

E^^\D,V)^ lim E,{Z,D,V) 

= lim [E{Z,D,V)-EMZ,D,V)] (3) 

Z—>(yo 

= E^^\D,V)-E'il{D,V). 



Despite intensive study PITO], the coefficient E^^^D, V) 
has not yet been reported in closed form. However, the 
accurate numerical estimate 



has been determined for the important D = 3 case [lOj . 
Combining (|4| with the exact result [8j 
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yields a value of 

e!:^^ = -0.046 663 253 999 48 



(5) 



(6) 



for the helium-like ions in 3-dimensional space. 

In the large- Z? limit, the quantum world reduces to 
a simpler semi-classical one |11| and problems that defy 
solution in D = 3 sometimes become exactly solvable. 
In favorable cases, such solutions provide useful insight 
into the D = 3 case and this strategy has been success- 
fully applied in many fields of physics [TH US]- Indeed, 
just as one learns something about interacting systems by 
studying non-interacting ones and introducing the inter- 
action perturbatively, one learns something about D — 3 
by studying the large- £* case and introducing dimension- 
reduction perturbatively. 

Singularity analysis [H] reveals that the energies of 
two-electron atoms possess first- and second-order poles 
at £) = 1, and that the Kato cusp [El [16] is directly re- 
sponsible for the second-order pole. In our previous work 
[T71 ITS] , we have expanded the correlation energy as a se- 
ries in 1/(D — 1) but, although this is formally correct if 
summed to infinite order, such expansions falsely imply 
higher-order poles at -D = 1. For this reason, we now fol- 
low Herschbach and Goodson [Tni H^ , and expand both 
the exact and HF energies as series in 1/D. Although 
various possibilities exist for this dimensional expansion 
[m [2TVE5] , it is convenient to write 



e^^\d,v) = 
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TABLE I. _B^2,o) ^ 
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coefficients for 



various systems and v{r) — 1. 
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Helium 


-1 


5/8 


1/2 


1/8 


0.424479 


Airium 
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7/24 


1/6 


1/8 


0.412767 


Hookium 
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1/4 


1/8 


1/8 


0.433594 


Quartium 
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5/24 
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1/8 


0.465028 


Sextium 
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3/16 
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0.486771 


Ballium 
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0.664063 



where 



ij(2,l)(^)=^(2,l)(^)_^(2,l)(V^) 



(10) 
(11) 



Such double expansions of the correlation energy were 
originally introduced for the helium-like ions, and have 
lead to accurate estimations of correlation [Ml |^ and 
atomic energies [3S1 ^7\ via interpolation and renormal- 
ization techniques. Equations (I?]), (|8| and ^ apply 
equally to the ^ S ground state of any two-electron system 
confined by a spherical potential V{r). 

For the helium-like ions, it is known fTU\ [5D1 HB] that 
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384' 
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and we have recently found J7] that Ec ' (V) takes the 
same value in hookium (two electrons in a parabolic well 
P5H5^ ). spherium (two electrons on a sphere [33H3S]) 
and ballium (two electrons in a ball [T51 1371 I3H])- In 
contrast, we found that Ec ' (V) is F-dependent. The 
fact that the term Ec is invariant, while -Ec ' varies 
with the confinement potential allowed us to explain why 
the high-density correlation energy of the previous two- 
electron systems are similar, but not identical, for D = 3 
[T71[TB]. On this basis, we conjectured [T7] that 



Ei'HD,V) 



1 
81)2 



C{V) 



(13) 



holds for any spherical confining potential, where the co- 
efficient C{y) varies slowly with V{r). 

In this Letter, we prove that Ec \s indeed universal, 
and that, in the large-D limit, the high-density correla- 
tion energy of the ^S* ground state of two electrons is 



given by ( 13 1 for any confining potential of the form 



\/(r) =sgn(m)r™u(r), (14) 

where v{r) possesses a Maclaurin series expansion 

v{r) = uq + vir + V2r'^ + (15) 



In order to prove the conjecture ( 13 ), we start with the 
conventional Schrodinger equation 



and the general Hamiltonian 



(16) 



H = -\{^\+Vl) 



zm+2 



[V{n) + V{r2)] + --, (17) 



ri2 



where Z is the confinement strength and ri2 = |ri — r2\ 
is the interelectronic distance. After the Jacobian- 
wcighted transformation 



^ = rf-irf -isin^-2, 



(18) 
(19) 



where is the interelectronic angle, the Schrodinger equa- 



tion (|16|) becomes 



Ed'^d, 



(20) 



in which, for states with zero total angular momentum, 
the kinetic, centrifugal, external and Coulomb operators 
are respectively 
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We now need to recast the Schrodinger equation so 
that perturbation theory can be applied. To achieve this, 
we successively introduce the scaled quantities 



A 

7 

kZ 



z 



where k = A™+2 , and introduce the scaled energy 
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The Schrodinger equation then takes the simple form 

jf + U + V+^w]<^D^£D'i>D, (28) 

and it is clear that perturbation theory can now be used 
to expand the energy both in terms of Z and A. 



In the D = oo limit, the kinetic term vanishes and 
classical electrostatics cause the electrons to settle into a 
fixed ( "Lewis" ) structure [19] that minimizes the efi^ective 
potential 



X =U + V 



Z' 



-w. 



The minimization conditions are 

dX{n,r2,0) _ dX{n,r2,i 



dri dr2 

dX{r,,r2,e) 



= 0, 



0, 



(29) 

(30) 
(31) 



and the stability condition implies m > —2. Assuming 
that the two electrons are equivalent [32], the resulting 
exact density and energy are 

l^ool' = S{r, ~ r^)Sir2 - r^)S{e ~ 9^), (32) 

£oo ^ X{roo,roo,Ooo), (33) 

where 5 is the Dirac delta function. Substituting Taylor 
expansions of Too and 
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1 



into (30) and (311 yields 
1 Vl ~ 



TO + 2 \2V2 
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where a"*^™"'"^-' = sgn{m)mvo. The m = case requires 

special attention, and is found by taking the ttt, — > limit. 

For the HF energy, things are simpler. The HF wave 

function is independent oi 9, so the only angular depen- 



dence comes from the Jacobian (19). Moreover, because 
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it follows [2Q] that 9^ = 7r/2. Solving (l30|, one finds 



that r^ and r^o are equal to second-order in 1/Z. Thus, 
in the large-D limit, the HF density and energy are 
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)<5(r2-r^^)<5(0--), (37) 
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and correlation effects originate entirely from the fact 
that 6*00 is slightly greater than 7r/2 for finite Z. 

Expanding ( 33 ) and ( 38 ) in terms of Z and D yields 



E^^-'>HV) = -- - 
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E- 



(2,0) 
HF 



iV) = - 



8 2(m + 2) ' 

1 
2(m-h2)' 



(39) 
(40) 



thus showing that both E^"^'^^ and E^^ ' depend on the 
leading power m of the external potential but not on v{r). 
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FIG. 1. Coefficients of the exact (dashed), HF (dotted) and 
correlation (solid) energies with respect to m, for v{r) = 1 
(Eqs. ([7|, ([8f and ^). 



Subtracting these energies yields 



(41) 



and completes the proof that, in the high-density limit, 
the leading coefficient _Ec of the large-D expansion of 
the correlation energy is universal, i.e. it does not depend 
on the external potential V{r). 

What is the origin of the constant in Eq. ([41])? It 
comes directly from the leading coefficient (l/4-\/2) in 
the 1/Z expansion of 9ryo (Eq. ([35])) and, because that is 



determined via Eq. (31 ), it is independent of the external 
potential V{r). This reveals that Eq. (41) applies to a 



pair of electrons in any radial external potential, but not 
to anisotropic external potentials. 

Detailed analysis of -Ec shows that it results from 
contributions of -t-1/8 and —1/4 from the centrifugal po- 
tential U and the Coulomb operator VV, respectively. The 
external potential V, which contributes identically in the 
exact and HF treatments, does not contribute to the cor- 
relation energy. Kato has made a similar argument |15j 
to explain the behavior of the wave function as ri2 — >■ 0. 
In a D-dimensional space, the Kato cusp condition is |16j 
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dri2 
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'■12=0 



D-1 



*D(ri2 = 0), 



(42) 



and arises from the cancelation of the singularities in the 
Coulomb operator and the D-dependent angular part of 
the kinetic operator [B]. These observations suggest a 



connection between the result (41 1 and the Kato cusp 



(42). For large but finite D, the discovery that the Kato 



cusp plays a key role in the large-Z limit would not be 
surprising for, in this limit, the only relevant information 
is the behavior (42) of the wave function near ?'i2 — 0. 



The E^^'^^ and E^^ ' coefficients can be found by con- 
sidering the Langmuir vibrations of the electrons around 



their equilibrium positions [ini [20] ■ The general expres- 
sions depend on vq and vi, but are not reported here. 
However, for v{r) = 1, which includes many of the most 
common external potentials, we find 
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85 
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9/32 



(to + 2)3/2 

1/2 , 



f/f6 



(to + 2)1/2 (m + 2)1/2 + 2' 



(43) 



showing that _Ec ' , unlike _Ec ' , is potential-dependent, 
ft is singular at m — —2, tends to —85/128 as r?T. ^- oo, 
and reaches a maximum of —0.388482 at to. « —0.344223. 
The latter value of to corresponds to the minimum of the 
correlation energy in the large-D limit. Numerical val- 
ues oi Ec are reported in Table |f| for various systems, 
and the components of the correlation energy are shown 
graphically in Fig. [T] 

In conclusion, we have proved that the leading term 
Ec ^ — 1/(8I?2) in the large- Z? expansion of the high- 
density correlation energy of an electron pair is invariant 
to the nature of the radial confining potential. Although 
formally divergent J40] , truncated l/D expansions have 
been found to be a powerful tool for the exploration of 
correlation effects and, in the present study, they help to 
explain the observation that, in finite-dimensional spaces 
such as D = 3, the correlation energy depends only 
weakly on the confining potential. 
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